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Abstract 

We combine, in a single set-up, the complex time parametrization in path inte- 
gration, and the closed time formalism of non-equilibrium field theories to produce a 
compact representation of the time evolution of the reduced density matrix. In this 
framework we introduce a cluster-type expansion that facilitates perturbative and non- 
petrurbative calculations in the realm of open quantum systems. The technical details 
of some very simple examples are discussed. 

PACS: 03.65.-w, 03.65.Yz 



1. Introduction. 



In recent years there has been increasing interest in the consistent description of the 
dynamics of open quantum systems [TfjH]. Quantum decoherence and dissipation are very 
important phenomena in many different areas of physics. A non-exhaustive list includes 
problems from quantum optics to many body and field-theoretical systems. Dissipative 
processes play a basic role in the quantum theory of lasers and photon detection, and they 
are equally important in nuclear fission and the deep inelastic collisions of heavy ions. More 
recently, the influence of the environment on a quantum system emerged as an issue of 
crucial importance, not only due to its fundamental implications, but also due to its practical 
applications in quantum information theory (8 10 . In fact, during the last decade, many new 
discoveries regarding the physics of open quantum systems were made. Primary examples 
of a promising progress can be found in the rapidly developing field of quantum optics and 



the connected continuous variable systems in quantum computation 12 , 13 



Theoretical studies of decoherence and dissipation in quantum mechanics are centered 
on the time evolution of the reduced density matrix of a system embedded in a specific 
environment. The basic tools for studying the reduced dynamics are either effective equa- 
tions of motion, where the dynamics of the environment are eliminated, such as the Lindblad 
master equation [6}|7), or the influence functional technique introduced by Feynman and Ver- 
non 



14 . The latter is based on the path integral approach, and was used by A. Caldeira and 



A. Leggett 15 in the study of the quantum Brownian motion more than twenty years ago. 
In most cases, however, neither the Lindblad equation nor the influence functional can be 
exactly evaluated, since the interaction between the system and the environment is too com- 
plicated. In fact, the simulation of the environment by a system whose degrees of freedom are 
treated as random variables following a more or less simple distribution, is a rather common 
practice. Therefore, one usually relies on some simple, specific system-environment models: 
a harmonic oscillator or a two-level quantum mechanical system embedded in a (thermal) 
bath of other harmonic oscillators or other spin systems. In the present work we aim to 
introduce and investigate calculational tools capable of exploring the behavior of an open 



2 



system in interaction with a specific quantum environment. To be precise, we investigate 
the possibility to extend the calculational capability of the Feynman- Vernon path integral 
approach by adopting and combining definite functional methodological tools already known 
from different research fields. The first such tool is a combination of the well-known "closed 
(real) time formalism" |18| with the (equally well-known) imaginary time formulation [2] in 
the context of path integration. The compound result, called "closed complex time formal- 
ism" (or CCT ), enables us to isolate, in a simple and compact expression, the influence of 
the environment on the evolution of the system. It is well known that, in general, the inte- 
gration of the environmental degrees of freedom does not produce a local "effective action" 
that controls the dynamics of the sub-system. The so-called Feynman- Vernon action, which 
incorporates the influence of the environment, is a highly non-local object: it is non-local 
in time and in space. The proposed CCT technique has a well-defined result: it produces 
an influence functional that can be viewed as an action local in space. In this action the 
paths are defined on the complex plane and they are parametrized with the help of a "time" 
running on a specific contour of the complex plane. The interest in such formalism is not 
"theoretical" but practical: one hopes to transfer the existent richness of perturbative and 
non-perturbative path integral techniques into the realm of open quantum systems. 

Our second suggestion, strongly related to the first one, is the application of the so-called 
"cluster expansion" in the CCT context. The foundation of the application of this very pow- 
erful technique, of course, lies in the spatial locality -on the complex plane- of the influence 
functional. The cluster -or cumulant- expansion results to an expression that can be viewed 
as the "effective action" that governs the dynamics of the system after the elimination of the 
environmental degrees of freedom. However, in general, the cluster expansion produces an 
infinite series that contains all the orders of the environmental connected correlators and, if 
it is to be useful, some kind of truncation is necessary. As a first step in this direction, we 
consider the case in which the environmental correlators are of very fast decrease. Our for- 
malism allows us to prove quite generally and without any reference to a specific model, that 
the two-point environmental correlator (which is the most important in our approximation 
scheme) has all the properties that can lead the subsystem to decoherence and dissipation. 
It is worth noting that our proposal can be extended to systems with an infinite number 
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of degrees of freedom, such as the electromagnetic field interacting with matter or other 
field-theoretical systems. 

The remainder of the paper is organized as follows: In Section 2 we present the details of 
the complex time formalism in the context of the path integral formulation of the Feynman- 
Vernon influence functional, and we discuss the assumptions under which the aforementioned 
formalism is applicable. In Section 3 we apply the cluster expansion in the framework of 
the CCT formalism, and we discuss the emergence of some quite general and very important 
properties of the influence functional. In subsection 3.1 we provide a specific example of an 
environment which is just a simple harmonic oscillator (or a collection of non-interacting 
harmonic oscillators). In Section 4 we consider the case of an environment in which the 
correlations decay very fast after some characteristic time interval. This stochastic behavior 
truncates the cluster series, enabling explicit calculations pertaining to the open system 
per se. As a first step in this direction, in the same section we calculate the entanglement 
entropy of a simple harmonic oscillator. Finally, in Appendix A we present the details of the 
calculation needed for deriving the results appearing in section 4. 

2. Time Evolution and the Closed Complex Time Formalism. 

The best way to interpret the usefulness of the closed complex time methodology (CCT 
from now on) is the examination of the time evolution of the reduced (environment averaged) 
density matrix of an open quantum central system (s from now on), which interacts linearly 
with its environment (e from now on). Adopting the usual starting point we assume that 
the total Hamiltonian can be written as the sum of two parts that refer to the system and 
the environment respectively, and a third part describing their interaction: 



The total system evolves in time unitarily and, consequently, the reduced density matrix 
changes in time according to the equation: 



The dynamical content of the last expression is incorporated into a time evolution operator 
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H = H s + H e + Hj 



(1) 




(2) 



that contains the degrees of freedom of the whole system: 



i 



t 

i uYj./\ 



tf(t,t ) = Texpj--jf dt/H(t')j (3) 

In the last expression we have taken into account a possible time dependence of the Hamil- 
tonian. Physically we understand such dependence in various ways; for example, we can 
imagine that, after a sudden quench, the coupling between the central system and its en- 
vironment changes to a different value, remaining constant henceforth. A case of physical 
interest arises when the coupling changes continuously and slowly enough to consider the 
evolution of the whole system as adiabatic. Another example is the well-studied case of an 
external time dependent field coupled linearly to the central system. In any case the operator 
T takes care of the needed time ordering. For now let us assume that, at the initial time to 
(for the sake of convenience, in what follows we shall assume that to — 0), the total system 
is prepared in a pure disentangled state 

IV> = hk)®hk). (4) 

Consequently, we can rewrite the reduced density matrix in the form: 

p R (t) = tr e [U(t)p s (0) <g> p e {U)U\t)\ (5) 

Denoting by x and q the coordinates of the central system and the environment respectively, 
and by X = [x\, Xd', qi, Qd) the coordinates of the whole system collectively, eq. ^ 
can be written in the well known form: 

p§ g (t) = J d D x" J d D x"' P l„ slll (0)J(x, g , x", x" f ; t) (6) 

where the propagating kernel can be read from the expression: 

J{x,x',x",x"';t) = J d D q J d D q" J dVV^W^glC)'^)!-^'')^'''!^^)^)- (7) 
Our next assumption is that the environment is initially in its ground state: 

\A) = |0e) (8) 

Then, it can easily be shown that [19] : 



l - J I exp | - i | ° drC™ [g (3) ] - ^ jH dr4*> [g (2) ] I (9) 



q('i)(-0)=q" g< 2 )(+0)=g 1 



{ F) 

In the last expression we denoted by C e the Euclidean version of the Lagrangian describing 
the dynamics of the environment. The origin of eq.(|9]) can be traced back to the propagator: 

G e (q',t';q",t) = ^(<f | exp <^ -Ut' -t)H e \ \n e )(n e \q") = 



(10) 



Introducing the Euclidean time r = it, taking the limits, r = — Te, t' = 0, Te — > oo and 
assuming that the ground state is unique one can easily deduce that 



G e {<?\U-e-T E ) = (f\0 e )(0 e \q"')e~ T ^ /h 



1 + ( e ~ Tls ( Ene ~ Eoe)/h ) 



:n) 



and, consequently the ground state wave function can be determined through an integration 
of the Euclidean propagator: 

-o 



(«f |0 e ) ~ J dq"'G e (q", 0; f, -oo) = J Vq exp |-i | dr^fg] 



(12) 



g(-0)=«f' 

The above relation is the basis of eq. ([9]) in which we also introduced the normalization factor 

Z e = jvq exp|-^ + °°dr£f)[g]} (13) 

<f(-oo)=(f(oo) 

ensuring that £r e [p e (0)] = 1 and we used a numbering convenient for our future considera- 
tions. To proceed further we write: 

£( 4 ) (t)=x' 

{ x-,q]U(t)\X") = 

xW(0)=x" 

and 

xW(t)=x 

(X"'\U(t)\X) = I PjeW^^-gl^Hoj-riexpl^^i'^Wlj. (15) 

x(i)(0)=x"' 



^'(<)=f f t ~) 

#(£)|X") = J VX^ S[q< 4 \t) - q\5[q< 4 \0) ~ exp U jf ^£ e [X( 4 )] I (14) 



{0,+iOh 



(1,+W) 



(0,-/0) 



Figure 1: The Contour C 



Inserting eqs.Q, (14) and (15) in expression we find: 

T I -» -»/ ^// -»///. ,\ 

rj I 3/ . 2/ . Ob ■ Ob . L J 



xW(t)=x x^(t)=x' . 

" f J" dt'C s [x {4) } + l I dt'C s [x w ] j>^[x (4) ,x (1) ;t]. (16) 

gW(0)=S"' xW(0)=x" 

The last factor in the above equation defines the well-known Feynman- Vernon functional 



ly^y^j-; 
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which incorporates the influence of the environment to the time evolution of the system: 

4 



x exp 



h 



-o 



h 



h 



x 



(17) 



Up to this point the only difference of the last result from the usual line of thinking |1}]4] 
is that we consider the environment not as a heat bath in thermal equilibrium but as a 
quantum system -probably a very complicated one- in its ground state. 

The expression for the influence functional can now be considerably simplified if we 
introduce the complex variable z defined on the contour C shown in Fig. 1. This contour 



consists of 4 different straight lines: The first line L\ goes parallel to the real axis from the 
point z = t — zO to point z = — iO. The second line L 2 begins from the point z = — iO and, 
following a path along the imaginary axis, goes to z = — zoo. The line L 3 traces a path 
along the imaginary axis and joins the points z = + zoo and z = + iO. The last part of 
the contour is the straight line L 4 : It goes parallel to the real axis from the point z = + iO 
to point z = t + iO. It is now easy to be proved that the "action" in the influence functional 



(17) can be written as follows: 

S= / dz{£ e [q c ] + Ci[q c ,x c ]} = / dz{C e [q c ] + g c q c x c } . (18) 
Jc Jc 

The notation in the last equation is defined as follows: Along the lines Li, i = 1, . . . ,4 , we 
have written {xl v qij = (x^\ q^ 1 ) and we have introduced a contour dependent coupling g c 



with g Ll = g Li = g and g^ 2 = g>£ 3 = 0. In expression ( 18 ) we have explicitly assumed that the 
interaction between the system and the environment is linear and has the minimal coupling 
Lj = gxq. In what follows we shall also assume that the coupling g is time independent, but 
our considerations can easily be generalized to a time dependent coupling. 



To confirm that eq.(18) does indeed represent the action in the influence functional, let 
us note that along the lines L\ and L4 we can write z= t' — iO, and consequently: 

/ dz{C e [q Ll ] + Ci[q Ll ,x Ll ]} = 
Jli 



and 



/ dz{L e [q Ll } + Cj[q Ll ,x Ll }} = dt'L e+I [X^] (19) 

Jt-iO Jt 

dz{C e [q L4 ] + £/[g L4 ,£L 4 ]} = 



>L 4 

pt+iO pO 

/ dz{L e [q LA ]+C J [q L „x LA ]}= dt'L e+I [X^} (20) 

JO+iO Jt 



'0+iO 

Along the lines L2 and L3 we write z = — ir and thus: 

-0— ioo ^00 



n pi) — ZOO ^00 

dzC e [q L2 )= dzC e [q L2 }=i drCf\q^) (21) 

JL 2 JO-iO JO 



and 

"0+i0 /•- 



p pu+tu p— u 

dzC e [q Ls ]= dzC e [q L ,] = i drCf\q^} (22) 

J L3 JO+ioo J — 00 
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Inserting eqs.(19), (20), (21) and (22) into eq.(17) and imposing the boundary condition 



(f c (-zoo) = g c (+ioo) (23) 
we get the following compact expression for the influence functional: 

JF\x^\x^^t] = — [ Vq c exp <| ^ / dzC e [q c ] + ^ / dzg c x c q c 



q c (t+iO)=q c (t-iQ) 



C 



= ^exp |^ J dzg c q c x c jj = exp |^SW[^c]j (24) 

(Spy stands for the Feynman- Vernon action). 

As it is obvious from the above expression, the introduction of the complex time z defined 
on the contour C, has enabled us to interpret the influence functional as an integral over 
continuous paths with periodic boundary conditions. In fact, the compactness of the result 



indicated in eq. (24) is the essence of the CCT formalism. At this point it may be useful to 
summarize our assumptions. The first one was that initially the system and its environment 
were disentangled (see eq.Q). This assumption is not crucial either for the appearance of 
the influence functional or for the implementation of the closed complex time formalism. 
The basic assumption for the latter is that the time evolution begins from a ground state 
(see eq.(j8])). To confirm this statement, let us assume that initially the central system and 
the environment were entangled and the whole was in the ground state of the Hamiltonian 
H(P,X,g(0)). Consider now the time evolution with a Hamiltonian H(P, X , g(t)) in which 
the coupling between the system and the environment changes very slowly (or has a different 
constant value). The evolution of the reduced density matrix reads now: 

PM) = [d D q [ d 2D X" [ d 2D X"'p^ fll (0)(x' ,0(t)\X")(jt m \fr(t)\jt). (25) 



Following the reasoning that led us to eq. (9) we can write the initial density matrix in the 
form: 

= \ J ^ (3) / M (2) ex P { - \ J~° drC^[X^,9(r))- 



X( 3 ) (-0)=X" XW {+0)=X'' 



jH 'drC^[X^,g(r)]\. (26) 
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Using now the expressions (14) and (15) we arrive at the result: 



p* £ (t) = iy VX& J VX® J VX™ J VX^5[X^(0)-X^ 3 \-0)]: 



(t)=S> 



x8[X {2 \+0) - X (1) (0)]5[g< 4) (t) - ^ l \t)} exp 



h 



dt'C[X^\ 



h 



-o 



drC^[X^}- 



drC^[X^} + ^- / dt'C[X^\ 
o nJt 

Using the CCT formalism the last equation can be rewritten: 

x c (t+i0)=x' . 

f$&) = \\ VX c {z)e^r- jdzC[X c {z),g c {z)\ 

x c (t—i0)=x 



(27) 



Z 



x c (t+i0)=x" 



J Vx c (z)exp <~J dzC s [x c (z),g c (z)} I / exp J - J dzg c q c x c 

x c (t-i0)=x ^ J \ K C 

Thus the evolution of the reduced density matrix assumes the compact form: 



(28) 



x c (t+i0)=x" 

P§x(t) = Y V *c{z) exp i ^ / dz£ s [x c (z),g c (z)} + ^S FV }. 



(29) 



x c (t—i0)=x 

A detailed analysis of the time evolution of a non-product initial state, under a time- 
dependent Hamiltonian, will be presented in a forthcoming study. For the time being we 
focus on the case of a disentangled initial state and a time-independent Hamiltonian. 

3. The Cluster Expansion. 

It is self-evident that any further calculational step strongly depends on the dynamical 
details of the environment, as well as on the specific form of the interaction between the 



latter and the system. In any case, the compact formulation indicated in eqs.(24) or (29) 



can be combined with all the existent calculational technologies to produce concrete results 
in the field of open quantum systems. In this framework it is very convenient to use a 
well-known and very powerful technique: The so-called cluster or cumulant expansion. This 
fundamental technique is widely used in a great variety of problems from statistical physics 
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to quantum field theories [20]. The methodology has been extensively used in areas such 
as the resummation of perturbative series and non-perturbative estimations, among others, 
and has proven to be a very successful tool. 

In our case, the cluster expansion theorem can be read from the relation: 



oxi > \ ^ / dzg c q c x c 



where 



oo / • \ n 
I 



n=l 



h 



(>X I } ^ ( T ) / dzn- I dzi6 c {z n ,...,zi)K {n \z n ,...,z 1 ) 



K {1 \zi) = (Ci[q c (z\)\) q = g c {z 1 )x c {z l ){q c (zi)) 



(30) 



jr< a Wi) = (^to)]£/[?e(*i)]> 9 - (CMzdMCitoM]) 

= gc{z 2 )g c (z 1 )x c (z 2 )x c (z 1 ) (q c (z 2 )q c (z 1 )) q - (q c {z 2 )) q {q c {zi)) q 



K^\z 3 ,z 2 , Zl ) = (diq^Cjiqciz^dlq^)]),- 
-^(Ci[q c (z 3 )])q(Ci[q c (z2)]Ci[q c (zi)]) q + perms j+ 
+ (C I [q c (z 3 )](C I [q c (z 2 )]) q (C I [q c (z 1 )]) q , 



.etc. 



In eq.(30) we have introduced the chain of path dependent step functions 



@c( z n, ••■) z l) — @c( z n ~ z n-l)---@c( z 2 ~ z l) , 



(31) 



(32) 



which takes care of the time ordering needed whenever the variables z\ are integrated along 
the same line. The path dependent step functions that appear in the above expression can 
be defined with the help of a proper parametrization of the contour z = z(a), o e [0, 1] 
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with z(0) = t — iO and z(l) = t + iO. Since the time (real or Euclidean) flow follows different 
directions along different lines, we have introduced the following definition 

{9(a — a'), when C = L 2 a . 
(33) 
6 (a' — a), when C = L 13 

When the variables z are integrated along different lines the step functions become identically 
1 or 0: For example, if z G L\ and z' G L4 we define Ql x \jlX z — z') = 1 because the time 
along the line L\ decreases, and this happens after its growth along the line L4. 



The validity of eq.(30) with the definition (33) can be readily proven by expanding the 
corresponding exponentials. The proof can also be easily extended to the case of non- 
commutating quadratic matrices with the help of a proper time ordering. Taking into ac- 
count the above conventions, any well-known result of the ordinary path integration can be 



transferred into the complex time framework as it is defined by the expressions (24), (30) 



and (31). 



From the preceding analysis we saw that the influence of the environment has been 
incorporated into the correlators 

( l hT- 1 A^_ li (z n ,...,z 1 ) = (? C)i „(z n )...g Cjil (z 1 )) g (34) 

that must be integrated along a closed contour C defined on the complex plane and consisting 
of 4 lines in a definite order determined by the defining expression ^ for the evolution of 
the density matrix. The time flow along the aforementioned contour is not causal, in the 
sense that its growth (along the line L 4 ) comes after its decrease (along the line Li), a fact 
being taken into account in the properly defined path dependent step functions. 



As it is evident from the definition of the path integral in eq.(17), and the fact that the 
couplings disappear along the imaginary axis, non trivial correlations can exist only along the 
lines L\ and L4 or among them. This is closely related to the fact that initially the central 
system and its environment were disentangled. However, as we have already seen, the CCT 
formalism can also be applied if the system and its environment were initially entangled. In 
such a case, non trivial correlations can exist among all of the lines of the contour C. 



At this point, we can highlight the properties of the fundamental functions (34) by 



discussing some of the properties of the two point correlator which is supposed to be invariant 
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under space rotations and time translations: 

A'g! hil (z 2 ,z 1 ) = S i2il Gf(z 2 -z 1 ). (35) 

A first observation is that it must have a non vanishing imaginary part due to the imaginary 
period over which it is defined. To be concrete, let us consider the propagation along the 
line L\. 

G { l\(t 2 - t x ) = G R (t 2 - *0 + zG 7 (t 2 - tO. (36) 
Along the line L4 the time flow is reversed, and consequently: 

Gg(t 2 - tO ee G R (tt - t 2 ) + iGfa - t 2 ) = Gg(t! - t 2 ). (37) 

At this point we can appeal to the hermiticity of the density matrix: The influence functional 
must remain the same if we interchange x^' and x^ while taking the complex conjugate. 
The last action reverses the time ordering along the contour C, and consequently the function 
A ( L ^ must be anti-hermitian: 

iG%(t 1 -t 2 )]* = -G%(t 2 -t 1 ). (38) 



Thus we immediately conclude that the real part of the propagator (36) is an odd function, 
while its imaginary part is an even function of time: 

G R (t 2 - tO = -G R {U - t 2 ), G T {t 2 - tO = Gr{h - t 2 ). (39) 

The exchange contributions can also be deduced with the same reasoning: Since, as we 
have discussed, the time along L\ is after the time along L 4 , the exchange from the line L 4 

fry\ 

to the line L\ is controlled by a function G L ^ uLi (t 2 — ti) in which t 2 < ti, while the exchange 
from the line L\ to the line L4 must be controlled by a function G L ^ uLi (t 2 — t\) in which 
t 2 > ti. Clearly the relation 

Ggu^ta " h) = -[G ( H Ll {h - t 2 )]* (40) 

must hold. The trace of the reduced density matrix must be equal to one, and, consequently, 
the Feynman- Vernon action must go to zero as — > x^\ This can happen only if the 
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(forward) propagation L 4 — > L\ exactly cancels the (forward) propagation along L 4 , and the 
(backward) propagation L x — » L 4 exactly cancels the (backward) propagation along L\. 

GguL!^ " *i) = G «(*i - <2) - - *a) ( 41 ) 

and 

GtW*2 - *i) = -<^(*2 - *i) - iG/(t 2 - (42) 

These arguments show clearly that, quite generally, the order g 2 contribution to the Feynman- 
Vernon action assumes the form: 

= dt 2 j\u[xU{t 2 ) -x«\h)]Gr(t 2 - ti)[x (1) (ti) -x (4) (ti)]+ 

+ £ jf dfc dt![z (1) (t 2 ) - x( 4 )(t 2 )]G fi (t 2 - ^[xWfa) + xW(tO]. (43) 

It is now readily evident that the Feynman- Vernon action considerably changes the dynamics 
of the central quantum system. Its fluctuating part, which is connected to the imaginary part 
of the line propagator, reduces coherence. It is customary [TJ|2] and convenient to re-express 
its real part, which is connected to the real part of the line propagator, with the help of an 
even function 7(£ 2 — h) = "f(ti — h) through the relation: 

d 

Gnih-h) = —^(h-h). (44) 
at 2 

The function 7 introduces in the Feynman- Vernon action a term which, on the classical level, 



can be understood as a damping or "friction" term. Feeding eq.(43) with expression (44) we 
immediately find that: 

\sf y = J^*2 jTcftl[s (1) (*2) -^ 4) (t2)]G 7 (t 2 -* (4) (*l)] + 

+ k 1 1 dt2 r^i^ (1) fe)-^(t 2 )]7(t 2 -ti)[i (1) (^i) + i {4) (ii)]- 



JO 



-i 7 (0) / dh (^(^ 



+ ^[x (1) (0) + x (4) (0)] jf dt x [x^\t x ) - x (4) (ti)]7(^i)- (45) 



'0 
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At this point we must emphasize on the following observation: The last equation, being the 
exact contribution of the second cumulant in the cluster expansion of the Feynman- Vernon 
action, is not an approximate one. Despite the fact that it formally reproduces a colored 
-noise simulation of an uncontrollable environment [lj[2], it is the first term in a systematic 
approximation of the environmental dynamics. 



3.1. A Simple Example. 



As a specific example, let us compute, in the framework of the preceding analysis, the 
influence functional for the case in which the environment is just a simple harmonic oscillator 



m e U2 1 



m e u p q 2 . 



(46) 



In this very simple case only one term appears in the rhs exponential in eq.(24): 



h 



S 



FV 



h 



dz 2 g c (z 2 ) I dz 1 g c (z 1 )x c (zi)x c (z 2 )9 c (z2- zx)gP(z 2 - zi). (47) 



>c Jc 

The Green function appearing in the last equation obeys periodic boundary conditions and 
assumes the well-known form 

1 cos[uj e (\z 2 - z\\ c -T [2)\ 



Gq\z 2 — Zi) 



2m P u P 



sin(a; e r/2) 



with 



\z 2 - z x \ c = (z 2 - zx)[d c {z 2 - zi) - 9 c (zi - z 2 )\. 
The period is obviously imaginary T = —2iTg, and consequently: 



(48 



(49) 



G c (z 2 — Z\) 



2m P u e 



[i cosoj e \z 2 — z\ L + sincug \z 2 — zA 



2m P uo P 



-i0J e \z2-Zl\c 



(50) 



Given that = gi 4 = g and gi 2 = gi 3 = we can split the integration in (47) as follows: 

S FV = g 2 (I 11 + I M + I u ). (51) 



Where we used the notation: 



In = gf I dt 2 J dtxOih - t 2 )x w {t 2 )x w {t x ) 



cos ui e (t 2 — ti) sino; e (t 2 — tx) 

l ; h 



2m P u P 



2m P u P 
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/ dt 2 [ dt 1 9{t 2 -t x )y{t 2 )y{t 1 )[iG I {t 2 -t 1 )+G R {t 2 -t l )\. 
Jo Jo 



(52) 



In the last integral we have connected the result pertaining to the specific choice (46) with 



the general result (36) through the relations: 

,2 



Gi(t 2 — 1\ 



g 



cosu e (t 2 - t x ), G R (t 2 -t x ] 



g- 



2m P u, 



sina; e (t2 — t x ). (53) 



The second term in eq.(47) reads: 



hi = g 2 / dt 2 / «(t2-ti)x (4) (t 2 )x (4) (ti) 
Jo Jo 



cos ui e {t 2 — ti) sina; e (t2 — t\ 
I h 



2m P u P 



2m e uj e 



[ dt 2 [ dt 1 e(t 2 -t 1 )x^\t 2 )x^\t l )[iG I (t 2 -t l )-G R (t 2 -t 1 )}. 
Jo Jo 



(54) 



With the same reasoning the last term in eq.(47) takes the form: 



iu = g 2 dt 2 dt x [e{t 2 -t x )^d{t x -t 2 )]x^\t 2 )x^\t x ) 



,cosu e (ti — t 2 ) sinci; e (t 1 — t 2 ) 



2m P u P 



2m P u P 



t ft 

dt 2 / dt x 6{t 2 - t0a: (4) (*2)a: (1) (ti)[iGj(t 2 - h) + G R (t 2 - t x )}- 
u Jo 



t rt 

dt 2 / dt x 9{t x - h)x^{t 2 )x m {t x )[iGi{t x - t 2 ) - G R {t\ - t 2 )\. 

ii Jo 



(55) 



Inserting eqs.(52), (54) and (55) into eq.(47) we confirm the general result (43) with the 



specific expressions (53) for the real and the imaginary part of the line propagator. These 



forms can be readily extended to the case of a collection of N harmonic oscillators: 



N 



g 2 Gj(t 2 - t x ) ->• V 9n cosu ne (t 2 - t x )\ 



u 1 



2m P u r , 



g 2 G R {t 2 -t l ) 



N 

£ 

n=l 



9n 



2m P uJr, 



sin[cu ne (t 2 



(56) 



The last expressions are obviously the T — > limit of the well known result for an envi- 
ronment which is a heat bath consisting of a collection of harmonic oscillators in thermal 
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equilibrium 
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4. The Stochastic Environment. 

The cluster expansion discussed in the previous section, helped us to interpret the 
Feynman- Vernon action, and consequently the influence functional, as an infinite series over 
all possible correlations among the environmental degrees of freedom. However, it is evident, 
that such an interpretation can be useful only if the infinite series can be truncated with 
negligible error. The case of weak coupling between the system and its environment is a first 
and obvious example; we shall not discuss this occurrence in the present paper, but is worth 
noting that the use of the cluster expansion facilitates the resummation of the (asymptotic) 
perturbative series. 

In the present study we adopt the hypothesis that the dynamics of the environment 
establish a characteristic time scale r e after which all internal correlations decay very fast: 

Gj(t) = Gj — ► 0, j{t) = 7 (-) — ► for \t\ > r e . (57) 



T P I \T, 



The scale r e appearing in our starting relations (57), is such a time interval that, when it 
elapses, the environment returns to its initial state. We shall also assume (lj[2j|8 10 11 that 
there exists a second distinct time scale t s , characterizing the interaction between the two 
parts of the entire system and, consequently, the evolution of the reduced density matrix, 
which is much larger than r e : r s ^> r e . 

In order to be more precise, let us assign an order of magnitude H-fT^H to the second 



order cumulant appearing in eq.(30). We shall consider as stochastic the limit: 



f 1 1^( 2 ) 1 1 -> 0. (5* 

lb 

As clearly shown by its definition ||10 2 )|| is a measure of the average "strength" of the 
interaction between the central system and its environment: \J ||iT( 2 )|| ~ (V). Defining the 



time scale r s as r s ~ h/(V), the limit indicated in eq.(58) can be obviously rephrased as 
tJt s -»■ 0. 

We can now examine how the cluster expansion is formed at the stochastic limit. Assum- 
ing that (q c )q = the first non- vanishing contribution comes from the second order term, 



which, following the discussion in the previous section, assumes the quite general form (45). 
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As we are interested for t ^> r e we take into account eqs.(57) and (58), and, performing 
the expansion 



we get 



g~ I dt 2 [x {1 \t 2 ) - x^\t 2 )}Gj{t 2 - t x ) « a[x^\t x ) - x^\t x )\. 



In the last expression we have introduced the quantity: 

/oo 
-oo 



(59) 



(60) 



In the same way, the second term in the rhs of eq.(45) can be approximated as follows: 

g 2 f dt 2 [x^\t 2 ) - x^\t 2 )Ut 2 - t x ) » \[xW(h) - x^ih)}, (61) 
Jo 



with 



\ = g 2 dt 2l {t 2 ). 



(62) 



With the help of a time rescaling ti = T e ti and using the defining relation for the 7 function 
(see eq. (43)) we can estimate that: 

A 



OC T> — > 0. 



a 



(63) 



After the preceding approximations the second order contribution to the Feynman- Vernon 
action reads: 

'0 2fc Jo 



« ^(2) _ O" 



h S FV = -^ J dt l [x^\t 1 )-x^{t 1 )f + ^ I dt l [x( 1 \t 1 )-xW(t 1 )][x ( - 1 \t 1 ) + xW(t 1 )}- 



-^7(0) J dti[(^(ti)) -(^ (4) (ti))' 



+ 



iA 

2n 



(x^(0)) 2 -(x( 4 )(0)) : 



(64) 



Our claim is that, at the stochastic limit (58), the cluster expansion, and consequently the 



Feynman- Vernon action, is dominated by the second order cumulant which, in this case, is 



expressed, by the above written eq.(64). Indeed, each of the terms in the cumulant 



expansion represents a cluster that must be integrated over time intervals much larger than 
the time scale characterizing its exponential decay. Thus in the integrals 



lM= / dt n / dt n ^... / dt l K^ n \t ni ... 1 t l ) 
Jo Jo Jo 



(65) 
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the main contribution comes from time intervals \t\ — ti\ ~ r e , i = 2,3, .... Expanding the 



integrand as we have done in eqs.(59) and (61), we conclude that: 

t£> = °©' < 66 > 

This conclusion can be used to give concrete meaning to the environment characterized as 
stochastic: It is the environment whose influence can be approximated by keeping only the 
second order correlator in the cluster expansion. 

In other words, the Feynman- Vernon action, at the stochastic limit, can be approximated 
as follows: 

S FV « SP V + o( T A. (67) 



At this point we must underline, once again, the strong resemblance of our result (64) 



to the case of the so-called Ohmic environment 1144 , 15 ; that is, to the case of the quantum 



mechanical simulation of a white-noise reservoir. Despite the fact that the expression (64) 
for the Feynman- Vernon action is, in both cases, formally the same, our result must be 
understood in a different context: It is the first term in a systematic approximation of an 
exact result which is supposed to be valid at zero temperature. The parameters appearing in 



eq.(64) are not phenomenological, but they are strictly related to the two-point correlation 
function of the environment, and, in principle, can be calculated at least numerically. In 



the same context, the expression (24) which is approximated by (67), does not represent 
the introduction of a random complex-valued Gaussian stochastic force: It is the specific 
environment under consideration and its dynamics that justify the stochastic approximation. 
Having in mind the extension of our work to infinite degrees of freedom, the non-Abelian 



gauge theories 21 constitute the primary example of such a stochastic behavior. 

In the present study, the undertaken task is, so to speak, "phenomenological" : given the 



approximation (67) for the influence of the environment, we try to estimate the consequences 
on the central system. 



In any case, the result (67) considerably facilitates the process of determining the time 
evolution of the reduced density matrix. The final result depends, of course, on the initial 
state of the central system, as well as on its specific dynamics. In what follows we shall 
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consider the case in which the central system begins from its ground state 



IV.) = |0 S >. (6* 

In such an occurrence we can use for p s (0) an expression analogous to the one (cf. eq.(|9 
used in the previous section for the environmental density matrix: 

pi ; ,.,„(0) = 1 /" Px( 3 ) / Pf( 2 ) exp I - i p° drCf >[x< 3 ) 

x(3)(-0)=x"x( 2 )(+0)=x"' ^ 
1 /-+oo 

x exp 

^ - J / dr/^^] j>. (69) 



Inserting the last expression into eq.([6]) we immediately get, at the stochasticity limit, the 
following path integral representation for the reduced density matrix: 

x c (t+iO)=x' , ,. 

PM) Vx c (z)exp \ l -JdzC s [x c {z)] + ^[x e (z)\ V (70) 

x c (t-iO)=x ^ ' 

As expressed in the last equation, the result for the reduced density matrix is simple and 
compact. This is due to the complex parametrization of the paths under integration. To 
obtain the final result, the integration over the central degrees of freedom must be performed 
and, obviously, this is a task that cannot be exactly accomplished in the general sense: some 



kind of approximation is needed. In any case eq.(70) sets the scene where any available 
approximation technique can be performed. We can demonstrate the calculational abilities 
of our formalism by considering,once again, the zero order approximation i.e., the simple case 
in which the system is just one simple harmonic oscillator (we neglect any space index): 

£s{xc{z)} = ^mx 2 c - ^mu 2 x 2 c . (71) 

It is now suffices to observe that the contribution from the Feynman- Vernon action is 
quadratic, and consequently, the dependence of the reduced density matrix on the boundary 
values x and x' can be deduced just from the classical path: 

/ J2 \ ro( 2 ) \ r cL] 

m \d? + U ) x ^ z) = 5x* (z) ' ** ( * + ' 0) = X '' X *' (t ~~ t0) = X - (72) 
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In the last equation the rhs must be read in terms of the stochastic limit (64). Thus we 
readily obtain: 



vm , 



f$ x (t) ~ exp { -^Wxf(t + zO) - xx? (t - zC 



X 



(73) 



The last two terms appearing in the rhs of the previous relation, cancel each other due to 
the quadratic nature of the truncated Feynman- Vernon action. Thus we conclude: 



P x x (t) ~ exp { ^[x'xf(t + zO) - xif (t - zO)] 



(74) 



The appearance of the classical trajectory in the last equation calls for the solution of the 
equation of motion (70). This is a lengthy but straightforward task, and it is presented 
in full detail in Appendix A. At this point it is enough to observe that the dependence of 
the classical solution on the boundary values x and x' is easily determined using the quite 
general ansartz: 



xf{t) = \u{t)x' + l -f3{t)x, x { 2(t) = ^(t)x' + l -5(t)x. 



(75) 



In the Appendix A we determine the coefficients in the above relations and confirm the 
validity of the relations S(t) = a*(t) and j(t) = (3*{t), which are necessary for the hermiticity 



of the reduced density matrix. Inserting expressions (75) in eq.(74) we find that: 



p * x (t) = C(t) exp { ^ x' 2 a{t) - x 2 a*(t) + x'x{(3{t) - /?*(*)) 



Ah L 



(76) 



The suppression of the off-diagonal terms in the representation (76) of the reduced density 



matrix is obviously related to the non-zero imaginary part of the function a(t), which in 
turn, as we confirm in the Appendix A, is related to the non- vanishing imaginary part of the 



environmental correlations. The normalization factor in equation (76) is now determined by 
demanding: 



C(t) / fix exp 



m 
2h 



x 2 %(a(t) +/3(t)) 



1. 



(77) 
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The explicit calculations presented in Appendix A show that 

S(a(t) + /?(£)) = (78) 

yielding the conclusion that C = 1/L — > 0, where L is the volume of the space in which the 
system lives. In this case the reduced density matrix reads: 

P*>M ~ exp j^a(t)(*' 2 - * 2 )} exp {^a(t)(x' - x) 2 ) . (79) 

The explicit form of the function a(t) is presented in Appendix A. Here suffice it to note that 
is a positive definite increasing function of time. It is strictly related to the imaginary 
part of the environmental second order correlator since oc a. Thus, the real factor of 



the density matrix (79) is formally the density matrix of a free particle in a heat bath of 
temperature fcgT = ^Qa oc a. 

The exact time dependence of the function a(t) is tied with the value of the quantity: 



A 2 7(0) 

q 2 = + 2^-u 2 . 80 

4m m 



If q 2 > 0, a(t) becomes time independent for t\q\ 3> 1 and 



9ta«— T^r, 3fe*w — + 2|g|. (81) 



m \q\ m 

For q 2 = and for {oj — \/m)t ^> 1, a(t) is again time independent: 

%a ^— 2 -—, $ta^2tu. (82) 

muj — X/m 

If q 2 = — k 2 < 0, and for kt ^> 1, remains an increasing function of time: 

^ , . a A; 2 + (w - A/2m) 2 . . 

w m [(u - A/2m)smH + k cos kt] 2 K ' 

The reduced density matrix is the crucial quantity for the physics of an open system, 
playing a key role for determining all the system properties. As an interesting example, we 
shall focus on the entanglement entropy 

S ent {t) = -tr s {p R (t)\np R (t)}. (84) 
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The calculation of the entropy can be performed with the help of the so-called replica 



method 19 . To apply it, one introduces the quantity 

tr s (p R r = f{n) = [ dx^ [ cfaW... f dx^p% xi2) p% x(3) ...p% xW . 



15) 



After calculating the function f(n) for integer n, we consider the function 

f{u) = tr s (p R r, u>0. 

Using analytic continuation we can find the entanglement entropy from the relation 

tr s {p R Y - 1 



(86) 



Sent. = - lim 

u-¥l V — 1 



-tr s {p H (t)\np R (t)} 



J7) 



Inserting eq.(79) into eq.(85) we get: 



(n „ \ n 

Hj dx«Jn«p{-^9fa(0( 



Consider now the propagation of a free particle with mass m from the point x to the point 
x' in the Euclidean time interval t% = 2/^sa(t): 



x(t E )=x' 

Vx exp 

x(0)=x 



m 

T 



dr x 2 (r) 



1/2 



exp 



m 
Ah 



%a(t)(x' 2 - x 2 ) 2 } • (89) 



Inserting the last expression into eq.(|88|) we find that: 

4ttH 



f(n) 



n/2 



nt t 



rn 



Vx exp { —— J dr x (t) 



(90) 



x(0)=a;(ni^) 

The last integral must be performed over periodic trajectories with period ntE- Thus we can 
immediately conclude that: 



f(n) 





n/2 


~m%a(t)L 2 ~ 


1/2 


m^sa{t)L 2 









(91) 



The entanglement entropy is now easily computed with the help of eq.(87): 

4nh 1 1 



Sent. — ~~ ^ ^ n 



m'<sa{t)L 2 



+ 



(92) 
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It is worth noting that, as it is well-known, the entanglement entropy S ent . ~ InL is not an 
extensive quantity: contrary to the thermal entropy, is not analogous to the volume of the 
space in which the subsystem lives. 

5. Conclusions and Perspectives. 

In this paper we have introduced two basic methodological tools for calculating the 
time evolution of the reduced density matrix and, consequently, the dynamics of an open 
quantum system. The first is the closed complex time (CCT) formalism, which combines 
two known approaches in a single set-up: The closed time formalism and the complex time 
one. This formalism enabled us to express the time dependence of the reduced density 
matrix, in terms of a compact path integral, in which the paths are parametrized on a closed 
contour on the complex plane. Our second metodological suggestion is the introduction of 
the cluster expansion which is a very powerful tool, tested in a variety of problems, where 
the environmental details can be successfully approximated by keeping only the two-point 
correlators. In this combined CCT-cluster expansion framework, we examined the case of the 
so-called stochastic environment in which the correlations are decaying very "fast" . In order 
to check our tools and examine the consequences of a stochastic environment, we performed 
a detailed "zero-order" calculation for the simple case in which the system is a harmonic 
oscillator. We found the explicit form of the reduced density matrix as a function of time 
and we calculated the entanglement entropy. Depending on the details of the environment, 
the entropy is either a constantly increasing function of time or an increasing function of 
time that saturates to a constant value. 

The purpose of this first work was to introduce and discuss the properties of a general 
formalism that can be applied in a variety of problems. We have confined ourselves only 
to a first -and in some sense trivial- application in order to demonstrate the underlying 
calculational machinery. In a forthcoming study we shall present the far more interesting 
case of the so-called quantum resonance. The general scene in such a problem is a double well 
embedded in a stochastic environment and in interaction with an external time dependent 
field. The path integral description of the tunneling and the role of the "classical" solutions 
in the framework of CCT is a very interesting and far from trivial problem that is under 
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investigation. 



Appendix A 

In this Appendix we shall determine the functions a(t) and (3(t) beginning from the classical 
equation of motion 

j2 \ ro(2) r d.l 

i , ,2 1 „cl./„\ _ UJ FVl x c 



Due to its nonlocal character the above equation must be examined independently in every 
segment of the contour C. 

Along the line L 4 the classical equation takes the form: 



rn 



where we defined 

mO 2 E-mu 2 + 27(0). (A.3) 
Along the lines L 3 and L 2 we have 

m(^-^)x^(r)=0 (A A) 



and 

m 



(^-- 2 )4?(-) = 0. (A.5) 
The last part of the classical equation refers to the line L x : 

m {^~ = ~ { X i - ia ) X ^ (A - 6) 

Seeking for continuous and differentiable solutions of the above system of classical equations, 
we impose the following boundary conditions: 

^{t)=x\ x«(0) = x®(0) 
x2 ) (-oo) = 0, x^0) = x<£{0) 
4 2 ?(0)=4 1) (0), *<?(+«)) = 

x$(t)=x, 4 ) (0)=4f ) (0) 
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and 



^?(0) = ^(0), i^(0) = -i^(0), iW(-oo) = iW(oo). 



,(2), 



(A.8) 



Equations (|A.4|) and (A. 5) can be readily solved with the help of the above indicated bound- 

4 3 V) = 4?(0)e^, 4? } (t) = x«(0)e— . (A.9) 



ary conditions: 



Using once again the boundary conditions (A.8), we find that: 



wx2?(0)=iS?(0) I c4 ) (0) = x«(0). 



Introducing the combinations 



,(±) 



the system of eqs.(A.2) and (A. 4) can be considerably simplified: 



d 2 A d 



(A.10) 



(A.ll) 



(A.12) 



^ mtf-^K^- - (A ' 13) 
The solutions j/ 1 ^ of the last equations are now trivially obtained and they lead us immedi- 
ately to the result: 



x 



(4) (t') = A lV f(t') + A 2 ^\t') + A 3 e a +*' + A 4 e- a - 1 ' 



(A. 14) 



x^(t') = A lV { l\t') + + A 3 e a+t ' + A A e- a - 1 '. 

In the above expression we have written: 

tpfW) = e a+t ' - 2i- [ dt"G(t',t")e a+t ", 
m Jo 



^\t') = e - a - t ' -2i- [ dt"G(t',t")e- c ' >' 
m Jo 



e " + * + 2i- f dt"G{t',t")e a+t ", 
™ Jo 



(A.15) 



(A.16) 



(A.17) 



(A.18) 
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^{ t ') = e- a - t ' + 2i- [ dt"G(t',t")e-° '" 
™ Jo 



(A.19) 



■±— + / A " , 2 ^(°) ,ji 
2m V 4m 2 m 

In eqs.(A.16) - (A.19) we used the Green's function 

- - tf) G(t f ,t") = -5(t' - t"), G(t,t") = G(0,f) = 
at u m at' J 

which assumes the form: 

e (a++a-)t/2-a + t" _ e -(a++a_)t/2+a_t" 



(A.20) 



(A.21) 



G(t', t") 



2(a + + a J) sinh[(a + + a_)t/2] 



e a+t ' - e- Q -*>(t" - t')+ 



+ 



e (a++a-)t/2-a-t' _ e -(a + +a-)t/2+a+f 

2(a_ + a+) sinh[(aq_ + a-)t/2] 



e~ a+t )6(t' -t"). (A.22) 



The coefficients in eqs.(A.14) and (A. 15) can straightforwardly be obtained with the help of 



the boundary conditions (A. 7) and (A. 10): 



A_(*) x' -x 
D(t) 2 ' 



(A.23) 



\+(t) x' — x 
D(t) 2 : 



(A.24) 



a_ + u x' + x //+(£) A_(£) — ji-{t)\ + {t) _ a t x' — x 



£>(t) 2 



(A.25) 



with 



A A (t) = 0+ - - X ' + - - / i +( t ) A -( t ) ~ V-MKjt) c a +t x' ~ X 



D(t) 2 



D(t)D(t) 



(A.26) 



£)(t) = A + (t)e- Q -' - A^(t)e a+ *, Z?(t) = (a+ - ^e'^ + (a_ + w)e 



(A.27) 



A ± (t) = #(0) + <p£'(0) - 2w, ^ = |(^(0) - v4 lj (0)). 



(A.28) 
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Inserting flA.22| and ( |A.23[ ) into ( |A.14[ ) and ( ]A.15[ ), we determine: 



;( 4 ), 



^+(t)A-(t)-M-(t)A + (t) e(a+ _ a _ )t 
5(f)I>(f) 



+ 



ct+(a:_ + w)e a+ * — — a;)e Q_ * 



Z>(f) 



(A.29) 



m = r± - _ (a+ _ a _ ) ^(t)A-(t)-/,-(t)A + ft) e(a+ - tt .) f+ 



£>(t) 



+ 



ct+(a:_ + w)e a+ * — a_(a: + — u)e a ~ l 



D(t) 



(A.30) 



7(*) 



A + (t)^ 1} (t)-A_(t)^(t) 
D(f) 



+ ( /. + (t)A_(t)-,_(t)A + (t) e(a+ _ a _ 



+ u)e a+t — — w)e" 



+ 



-a-t 



D(t) 

\ + (t)^\t)-\-(t)$\t) 
6(t) = ^ (a+ -a. 



(A.31) 



(a_)_— a_ jt I 



+ 



a+(a:_ + w)e a+ * — — a;)e Q_ * 



L>(t) 



(A.32) 



(The argument in all the functions is the instant t.) 

At this point we are ready to confirm some of the claims presented in the main text. We 
must distinguish two cases. The first is when: 



A 2 



> u 2 - 2 



7(0) 



4m 2 m 



(A.33) 



In such a case a± are real and consequently (p± = (<p± )*. Observing that A± = A±, 
fi± = — [i* ± we immediately see that: 



7 * = /?, 5* = a 



(A.34) 



and 



5a (t) = 



(A.35) 
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When 



A 2 2 7 (0) 
< u 2 -2-^ 



4m 2 



??? 



(A.36) 



we observe that a + = —a*_, (p± = (tp^)*, and since A±,/x± turn out to be the same as in 



the case (A.33), we verify once again the relations (A. 34) and (A. 35) 



When a± are real we straightforwardly obtain: 



3a(t) 



a 



m 



(a + -a-)t/"2 (a + -a-)t 

Ht) + -F»7*-Mt) 



D(t) 



D 2 (t) 



(A.37) 



3Ja(t) = 2^-]nD(t), 

(JjL 



with 



h(t) 



sinh[(a + + a-)t/2] 



1-e 



-(»++a_)t 



(a + + a_) 



+ 



(A.38) 



(a + + a_) 



and 



sinh[(a + + a_)t/2] 



te (a++«_)t/2 _ 2 Sinh[(a + + Q!-)t/2] 



(A.39) 



+ 



+ (a+ - w) 



sinh[(a + + a_)i/2] 



- t e -( a ++"-)</2 



(A.40) 



The last relations confirm that > 0. For ta± 3> 1 it is easy to check that ^sa and 
become constants: 



2a 1 

3a ~ , 9fta ~ 2a + . 

m at + + a_ 

The last relation holds as long as a + + a_ 7^ 0. If a + + a_ = that is if 

A 2 7 (0) , 
+ 2^^ = u 2 



we immediately find that 



a 



Am 2 m 



2t 



a 2 



m 1 + (u — A/ 2m) t^oo mco — X/m' 



(A.41) 



(A.42) 



(A.43) 



jto = 2 ^ + ^-A/2m)A/2m ^ ^ 
l + t(w-A/2m) *->oo 



(A.44) 



29 



When a± are complex we find that: 



+ 2(w - A/2m) sin 2 fe< + [k 2 - (u - A/2m) 2 ] I . (A.45) 

Using the fact x/ sin a; > 1 once again we can verify that 9fa(t) > 0. It also straightforward 
to see that: 

d — \/2m) — k 2 ] sin kt + k cos kt 

&a(t) = 2— InDft) = 2-^7 ^— A- , A.46 

w dt w (u - X/2m)smkt + k cos kt v ' 



where we have noted 

p_, 2 _27(0)_ 

m Am? 



;2 2 2 7(0) A 2 

k 2 =u 2 — -. (A.47) 
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